Magnetic penetration depth in single crystals of SrPd 2 Ge 2 superconductor 
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The in-plane magnetic penetration depth, A m (T), was measured in single crystals of SrPd2Ge2 
superconductor in a dilution refrigerator down to T = 60 mK and in magnetic fields up to Hdc = 1 
T by using a tunnel diode resonator. The London penetration depth, A, saturates exponentially 
approaching T — > indicating fully gapped superconductivity. The thermodynamic Rutgers formula 
was used to estimate A(0) = 426 nm which was used to calculate the superfluid density, p s (T) = 
A (0)/A (T). Analysis of p s (T) in the full temperature range shows that it is best described by a 
single - gap behavior, perhaps with somewhat stronger coupling. In a magnetic field, the measured 
penetration depth is given by the Campbell penetration depth which was used to calculate the 
theoretical critical current density j c . For H < 0.45 T, the strongest pinning is achieved not 
at the lowest, but at some intermediate temperature, probably due to matching effect between 
temperature - dependent coherence length and relevant pinning lengthscale. Finally, we find a 
compelling evidence for surface superconductivity. Combining all measurements, the entire H-T 
phase diagram of SrPd2Ge2 is constructed with an estimated H C 2(0) = 0.4817 T. 

PACS numbers: 74.25.-q, 74.25.Ha, 74.25. Sv, 74.25.Wx 



I. INTRODUCTION 

Superconductivity in the tetragonal ThCr2Si2-type 
SrPd 2 Ge2 was discovered first in polycrystalline 1 and 
later in single crystals 2 with the superconducting phase 
transition temperature (T c ) at 3.0 K and 2.7 K, respec- 
tively. The upper critical field {H c i) was estimated to 
be 4920 Oe at T = by using Helfand-Werthamer (HW) 
theory 3 based on the experimental data obtained only 
down to T = 0.7T C . 2 It has been found that T c and 
H C 2 can be slightly increased by chemical doping. 4 The 
London penetration depth and coherence length are re- 
ported to be A(0) = 566 nm and £(0) = 21 nm 5 and 
A(0) = 345 ± 30 nm £(0) = 25.6± 0.5 nm. 6 These 
values give the Ginzburg - Landau parameter of k = 
27 5 and k = 13. 5 6 , which makes SrPd 2 Ge2 a strong 
type-II superconductor. Furthermore, thermodynamic 2 
and tunneling spectroscopy measurements are consis- 
tent with a slightly strong - coupling s-wave Bardeen- 
Cooper-Schrieffer (BCS) superconductor with the zero- 
temperature value of the superconducting gap of Ao s=s 
2fcsT c 5 ' 6 , - not far from the weak coupling value of 1.76. 

This superconductor is interesting particularly be- 
cause of compositional similarity to the newly discovered 
isostructral Fe-and Ni-pnictide superconductors with 
comparable T c such as KFe2As2, BaNi2As2, and SrNi2P2- 
Although there is strong experimental evidence for nodal 
superconductivity in KFe2As2 (Ref. 7 and 8), the Ni- 
based ones have been shown to be fully gapped by ther- 
modynamic and thermal transport measurements. 9,10 
This naturally prompts the question: what the structure 
of a superconducting gap is in SrPd 2 Gc 2 ? So far, not 
much work has been done on SrPd 2 Gc 2 in this direction. 
Tunneling spectroscopy between 0.17T C and T c is con- 
sistent with a single, isotropic gap superconductor. 5 To 



verify this, however, the thermodynamic, thermal trans- 
port, and penetration depth measurements down to much 
lower temperatures are necessary to provide objective 
conclusions regarding the gap-symmetry in SrPd 2 Ge 2 . 

The magnetic penetration depth in superconducting 
state is among most useful probes to explore the super- 
conducting state. 11 In zero external magnetic field, it is 
defined by the London penetration depth. If measured 
with sufficient accuracy and low enough temperatures, it 
can be used to understand the angular variation of the 
superconducting gap on the Fermi surface. 11 ' 12 In the 
presence of vortices, the magnetic penetration depth is 
also influenced by the Campbell penetration depth which 
depends on the elastic properties of vortex lattice and is 
linked directly to the critical current density. 13 There has 
been only limited work performed in the mixed state of 
SrPd2Gc 2 6 ' 14 and no studies of the critical current den- 
sity over the full temperature and field range. 

Here we report precision tunnel diode resonator mea- 
surements of the magnetic penetration depth, A m (T), 
in single crystals of SrPd 2 Ge 2 (T c = 2.7 K) performed 
in a dilution refrigerator with temperatures down to 
T w 0.02T C and in magnetic fields up to 1 T w 2# c2 (0). 
The low - temperature variation of the London penetra- 
tion depth, AA(T), clearly indicates exponential satu- 
ration, and the analysis of the full - temperature range 
superfluid density is consistent with a fully gapped su- 
perconductor. The thermodynamic Rutgers formula was 
used to estimate A(0) = 426 nm which was used to calcu- 
late the superfluid density, p s (T) = A 2 (0)/A 2 (T). Anal- 
ysis of p s (T) in the full temperature range shows that 
it is best described by a single - gap behavior, perhaps 
with somewhat stronger coupling. The upper critical 
field, H c2 (0) = 0.4817 T was determined by the HW 
theory from the field - sweeps at different temperatures 
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down to 0.02T C . In finite magnetic fields, H < H C 2(0), 
the Campbell penetration depth, Xc(T,H), shows min- 
imum at an intermediate temperatures (rather than at 
the lowest temperature) which indicates non-monotonic 
variation of the theoretical critical current (j c ) calcu- 
lated from \c(T,H). Additional diamagnetic response 
detected above H c2 (T) is consistent with surface super- 
conductivity with H c3 (T) = 1.695iJ c2 (T). 15 



II. EXPERIMENTAL 

Single crystals of SrPd 2 Ge 2 were grown using a self- 
flux method as described in Ref. 2. The magnetic pen- 
etration depth was measured in a dilution refrigerator 
by using a tunnel diode resonator (TDR) technique (for 
review, see Ref. 11). The sample with dimensions 
(0.51x0.70x0.04) mm 3 with the shortest direction be- 
ing along the c— axis was mounted on a sapphire rod 
and inserted into a 2 mm inner diameter copper coil 
that produces rf excitation field with empty-resonator 
frequency of 17 MHz with amplitude H ac ~ 20 mOe, 
much smaller than H c \ of typical conventional supercon- 
ductors. Measurements of the in-plane magnetic pene- 
tration depth were done with both ti& c and H ac | c- 
axis. The shift of the resonant frequency (in cgs units), 
A/(T) = -G4ttx(T), where x(T) is the differential mag- 
netic susceptibility, G = foV s /2V c (l — N) is a constant, 
N is the demagnetization factor, V s is the sample vol- 
ume and V c is the coil volume. The constant G was 
determined from the full frequency change by physically 
pulling the sample out of the coil. With the character- 
istic sample size, R, Airx — (A/i?) tanh(i?/A) — 1, from 
which AA can be obtained. 11 ' 16 



III. RESULTS AND DISCUSSION 

A. London penetration depth 

Figure 1 shows temperature variation of the in-plane 
London penetration depth, AA(T), measured in a sin- 
gle crystal of SrPd2Ge2 superconductor which exhibits a 
very sharp superconducting phase transition at T c = 2.7 
K as shown in the inset, indicating a high quality, ho- 
mogeneous sample. In the main panel, AA(T) is shown 
with temperatures up to about 0.67T C . The saturation 
in T — > limit and almost flat temperature depen- 
dence, AA(T c /3) < 10 ran, indicate fully gapped su- 
perconductivity. Experimental AA(T) can best fit to a 
power-law function, AA(T) = AT n , with the exponent of 
n = 2.7 ± 0.1 and pre-factor of A = 12.2 ± 0.4 nm/K 2 7 . 
The fitting curve is shown in red solid line. A power- 
law function with such a high exponent has very weak 
variation at low temperatures, indistinguishable from the 
exponential behavior which is predicted for a fully open 
superconducting gap. In fact, the BCS low - tempera- 
ture form, AA(T) = A(0) y/nA /2k B T exp (-A /k B T), 
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FIG. 1. In-plane London penetration depth in a single crystal 
of SrPd2Ge2- Main panel: Open circles represent experimen- 
tal data. Solid and dashed lines represent power-law and BCS 
(single gap s-wave) low - temperature fitting. Dotted shows 
the data for KFe2As2 taken from Ref. 7 for comparison. In- 
set: London penetration depth in the full temperature range 
demonstrating a sharp transition at T c = 2.7 K 



where Ao is the maximum gap value at T = 0, fits the 
data equally well for T < T c /3 where it is expected 
to be valid. However, the best fitting is achieved with 
A(0) = 50 nm and A = 0.74k B T c . The latter is impos- 
sible in the single - gap clean limit where Ao ~ 1.76k B T c 
is expected. The value of A(0) is also much smaller than 
the reported value of 566 nm. 5 Similar low - temperatures 
features can be seen in two - band superconductors such 
asMgB 2 (Ref. 17), 2H-NbSe 2 (Ref. 18), Lu 2 Fe 3 Si 5 (Ref. 
19), and more recently LiFeAs (Ref. 20). However, as we 
show below, analysis of the superfluid density in the full 
temperature range is inconsistent with a two - gap clean 
behavior. Instead, it is more likely that we are dealing 
with moderate pair - breaking scattering (maybe due to 
well-known magnetic impurities in Pd) which results in a 
finite density of states inside the gap. We also point out 
that these temperature variation is very small compared 
to a known nodal superconductor KFe 2 As 2 (Ref. 7 and 8) 
with similar T c , but exhibiting much stronger tempera- 
ture dependence of AA, indicating significant amount of 
quasiparticles generated at the low temperatures, most 
likely due to nodes in the gap. 

For a metallic sample, the measured penetration depth 
above T c is determined either by the skin depth 5 or 
sample size. In case of skin depth limiting, the value 
of A(T > T c ) shown in the inset in Fig. 1 is one half 
of the actual skin depth. 21 Therefore, we can estimate 
normal state resistivity from the measurements using 
p = (2ttu;/c 2 )<5 2 . 20 For SrPd 2 Ge 2 with cj/2tt = / = 17 
MHz and 8/2 « 20 /jm, the calculated resistivity is ap- 
proximately 12 fiVL cm which is much less than the ex- 
perimental value of 68 jAl cm (Ref. 2). Therefore, we 
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conclude that it is in a sample - size limited regime. Us- 
ing the same equation, the estimated skin depth is 180 
nm which is close to half width of the sample. 

Finally, we note that the data exhibit a smooth transi- 
tion from superconducting penetration depth to the nor- 
mal state between T = T c and T* w 3.0 K, which has 
also been seen in transport measurement. citeSung2011 
Interestingly, T* = 3 K is the onset of superconduc- 
tivity observed in a polycrystalline sample. 1 Similar fea- 
ture has also been observed in a related superconductor 
BaNi 2 As 2 . 22 Perhaps this feature requires further study. 



B. Superfluid density 
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FIG. 2. Main panel: Calculated superfluid density, p s , with 
various A(0)'s. The slope of dashed lines was determined with 
AC P and \dH C 2/dT\T c using the Rutgers formula as described 
in the text. Inset: Variation of \dp s / dt\r c / A 2 (0) with varying 
A(0). Here p' a was determined from the experimental data. 
The gray band is a theoretical estimate with 5% hypothetical 
error in AC P and \dH C 2/dT\T c - 



While the low - temperature behavior is important, the 
superconducting gap can be probed at all energy scales by 
the analysis the superfluid density, p s (T) — A 2 (0)/A 2 (T), 
in the entire temperature range. 12 However, superfluid 
density requires knowledge of the absolute value of A(0). 
For SrPd2Ge2 A(0) = 566 nm was estimated using a 
dirty limit 5 and A(0) w 390 nm was extracted from field 
- dependent magnetization 2 within the London approxi- 
mation and A(0) ~ 345 ± 5 nm was estimated from the 
measurements of the field of first penetration. 6 So, the 
variation of the literature values is quite significant and 
we have to resort to another, thermodynamic, approach 
based on the Rutgers formula. In the Ginzburg-Landau 
regime, i.e. near T c , it can be shown that 



dp s 
dt 



Tc 



16tt 2 A 2 (0) 
b \dH c2 /dT\ Tc p 



(1) 



where 4>o = 2.07 x 10 G cm is a flux quantum and 
\dH C 2/ 'dT\x c = 0.26 T/K is determined experimentally 
(see Fig. 8). Specific heat jump AC P = 7381 erg/cm 3 K is 
taken from Ref. 4. Applying these thermodynamic values 
suggests \dp/dt\ T JX 2 (0) = 21.7 pmT 2 where t = T/T c 
is the reduced temperature. This quantity can be com- 
pared with the actual slope of calculated p s (t) with vari- 
ous A(0) at T c as shown in Fig. 2. In the main panel, the 
open symbols represent the calculated superfluid density 
with 300, 400, and 700 nm in triangle, circle, and square, 
respectively. The dashed lines are determined with the 
slope calculated by Eq. 1 for three values of A(0) quoted 
above. The line with 400 nm shows very good agreement 
with calculated p a while the line with 300 nm significantly 
underestimates, and the one with 700 nm overestimates. 
This procedure can be repeated with various values of 
A(0). The result is summarized in the inset where the 
solid triangle represents experimental slopes obtained by 
fitting experimental data near T c to a linear line. The 
gray horizontal band represents the theoretical value of 
\dp/dt\ T J\ 2 {0) = 21.7±2.2 pm~ 2 determined with a 5% 
hypothetical error in \dH C 2/dT\T c and AC p . In this way, 
A(0) can be determined at the intersection of the theo- 
retical line and experimental results, which provides that 
A(0) = 426 ± 60 nm that lies between the literature val- 
ues. With this value, the slope of p s at T c is determined 
to be -3.9. 
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FIG. 3. Calculated superfluid density, p s (T) = A 2 (0)/A 2 (T) 
using A(0) = 426 nm. Open circles represent the experimental 
data. The dashed dots and dashed lines represent single - gap 
weak - coupling s-wave BCS superconductor in clean and dirty 
limit, respectively. 

The calculated superfluid density with A(0) = 426 nm 
is shown in Fig. 3. The dot - dashed and dashed lines 



4 



show expectation for clean and dirty limit of a single- 
gap BCS superconductor in the weak - coupling limit, 
respectively. An attempt to use a two - gap (clean) 7— 
model 23 in the full - temperature range converges to a sin- 
gle - gap limit with A(0)/fcsT c =2.2. Therefore, the gap 
symmetry of SrPd2Ge2 is most likely represented by a 
single gap s-wave, perhaps with somewhat enhanced cou- 
pling strength. It was noted previously that the shape 
of p s (T) is rather close to a nonlocal - limiting case, ex- 
pected in type I superconductors such as aluminum and 
cadmium. 24 Similar argument was made in the work by 
T. K. Kim et al. in which SrPd2Ge2 appeared to be 
type I according to the intrinsic electronic structure de- 
spite the fact that experimental £(0) and A(0) put it in 
a strong type-II regime. 5 In any case, our study confirms 
that simple analysis with an isotropic Fermi surface is not 
sufficient and, perhaps, the results could be explained by 
taking into account a realistic band structure. We can, 
however, conclude that no nodes are present in the su- 
perconducting gap. 



C. Campbell penetration depth 

Figure 4(a) shows magnetic penetration depth, 
A\ m (T, H), as a function of temperature and magnetic 
field measured after cooling without magnetic field to 
target low temperature and then applying a DC mag- 
netic field of indicated amplitude (ZFC - solid lines) and 
upon cooling in field (FC - dashed lines). Increasing 
DC magnetic field not only suppresses the superconduct- 
ing phase transition and diamagnetic shielding, but also 
induces another diamagnetic phase between the normal 
state and apparent bulk superconductivity. This feature 
appears much clearer above H = 0.2 T and persists at 
least up to H = 0.6 T which is far greater than the 
upper critical field of the bulk superconductivity. This 
can hardly be understood with sample inhomogeneity or 
second phases since they both should affect the measure- 
ments in zero field. Systematic tracking of this feature 
reveals possible connection to surface superconductivity, 
which will be discussed later together with the general 
H-T phase diagram. In addition to temperature sweeps 
at fixed DC magnetic fields, magnetic field sweeps at dif- 
ferent fixed temperatures were performed, and supercon- 
ducting transitions are clearly detected as shown in Fig- 
ure 4(b). The determined bulk H C 2(T) is consistent with 
T C (H) determined from the temperature scans. We note 
that the anomaly above H C 2(T) was not detected in the 
field sweeping measurements. 

Another unusual experimental observation is ap- 
parently non - monotonic temperature variation of 
AA m (T, H) for 0.02 T < H < 0.4 T in both ZFC and 
FC data. The minima of AA m (T, H) are marked by open 
circles for clarity in Fig. 4(a), and the locations of these 
minima are nearly the same for both ZFC and FC data. 
This feature is too shallow to be due to paramagnetic im- 
purities (and should be most pronounced at H =0). 25 We 





FIG. 4. (a) Zero-field cooled (ZFC) and field cooled (FC) in- 
plane magnetic penetration depth measured in a single crystal 
of SrPd2Ge2 shown by solid and dashed lines, respectively. 
Open circles mark the minimum in A m (T) in each magnetic 
field, (b) Magnetic field sweeps at fixed temperatures. 



suggest that this behavior is cause by the size - match- 
ing of the temperature - dependent coherence length and 
relevant pinning centers. Furthermore, we observed the 
inversion of diamagnetic screening between ZFC and FC 
runs at about H — 0.2 T. The only way to obtain a differ- 
ence between FC and ZFC experiments in AC response is 
to assume a non - parabolic pinning potential 26 and our 
results would indicate that the effective pinning potential 
shape is strongly field and temperature dependent. Con- 
ventional superconductors, Al (Ref. 27) or YbSb2 (Ref. 
28) do not exhibit such upturns in a magnetic field. 

In the approximation of a linear elastic response of a 
vortex lattice to a small - amplitude AC perturbation 
(which is the case here), the total magnetic penetra- 
tion depth can be expressed as a sum of Meissner and 
vortex contributions, which are represented by the Lon- 
don and the Campbell penetration depths, respectively, 
= A 2 + A^. 13,29,30 Since we know A from the mea- 
surements in zero field, we can readily calculate Ac ■ The 
calculated Campbell penetration depth as a function of 
magnetic field shown in Fig. 5 appears sub-linear as ex- 
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FIG. 5. Isothermal Campbell penetration depth as function 
of a DC magnetic field calculated from the data of Fig. 4(a). 
Inset: AX C (T,H) = A c (0.06 K,H) - Xc(T,H) is displayed 
to show non-monotonic behavior. 



pected in conventional superconductors (~ v H) at a first 
glance. However, non-monotonic behavior is revealed 
upon closer inspection as shown in the inset. This non- 
monotonic behavior at low fields originates from non- 
monotonic temperature dependence discussed above. 




D. Theoretical critical current 

We distinguish theoretical critical current, which is a 
parameter entering the simple expression for the Camp- 
bell penetration depth, — j c = r p (f)o/\ c from the actual 
critical current that is affected by intervortex interac- 
tions and from the measured critical (persistent) current 
that is further affected by magnetic relaxation. Here 
r p is a characteristic radius of the pinning potential. 
We assume the r p to be equal to the coherence length 
£(T) = v/ 0o/27ri? C 2(T). Obtained theoretical critical 
current density is shown in Fig. 6, (a) as function of 
temperature at different fields and (b) - as a function 
of an applied field at different temperatures. The criti- 
cal current is non - monotonic as a function of temper- 
ature showing maximum at the intermediate tempera- 
tures. We attribute this to the matching between tem- 
perature - dependent coherence length and pinning land- 
scape. This assertion is plausible given monotonic behav- 
ior of the critical current versus field where the variation 
of the coherence length is much weaker. 31 Comparison 
with the previous works also re-affirms that we are deal- 
ing with the upper theoretical estimate of the critical 
current, approximately four times larger than estimated 
from the magnetic measurements. 14 On the other hand, 
direct comparison with our data shown in Fig. 8 shows 
that the maximum current density line found in this work 
is not an extension of the irreversibility line found in 



FIG. 6. Calculated theoretical critical current (a) as a func- 
tion of temperature in different magnetic fields and (b) as a 
function of magnetic field at fixed temperatures. Inset in (b) 
displays Aj c (T,H) = j c {T,H) - j o (0.06 K, H) to show the 
maximum j c being at a intermediate temperature. 



Ref. 14 and probably represents a crossover in the pin- 
ning mechanism reflected in the change of the effective 
pinning potential. 

Unfortunately literature data are only available for 
higher temperatures above T = 2 K and weak fields less 
than H = 0.1 T, so we do not know whether the non- 
monotonic behavior of j c (T) propagates to the relaxed 
persistent current density. 

Furthermore, our evaluation of the coherence length, 
£(T) (right axis in Fig. 7), allows experimental deter- 
mination of the temperature - dependent Ginzburg - 
Landau parameter, k(T) = A(T)/£(T), shown in Fig. 7 
(left axis). The result is compared with the Gor'kov 
theory 32 where temperature correction is introduced as 
A G (T) = k{T)/k(T c ). The trend is correct, but the mag- 
nitude of the variation is smaller then the observed. If we 
attempt to rescale the Gor'kov's result, the best agree- 
ment at the intermediate temperatures is achieved with 
k(T) = k (T c ) {2.6A G (T) - 1.6). Perhaps a better agree- 
ment could be achieved with a realistic band structure. 
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FIG. 7. Calculated Ginzburg - Landau parameter, n(T) = 
\(T)/£(T) (symbols, left axis) compared to the Gor'kov pre- 
diction, k (T) = k (T c ) Ag (T) (solid line) and to modified to 
best fit the data expression, k (T) = k (T c ) (2.6A g (T) - 1.6) 
(dashed line). Right axis shows estimated coherence length. 



E. H — T phase diagram 

Finally, a H-T phase diagram is established based on 
the measured AX m (T, H) as shown in Fig. 8. Open and 
solid circles represent bulk superconducting phase transi- 
tions determined by temperature- and field- sweep mea- 
surements, respectively. The red solid line throughout 
these symbols is the HW - fit 3 with H c2 (0) = 0.4817 T. 
This value was determined with the experimental data 
down to 0.02T C in our study, and it is very close to the 
estimated value H c2 (0) = 0.4920 T in Ref. 2 only with the 
the slope at T c by using the same theory. This is a good 
indication that H C 2(T) of SrPd2Ge2 is orbit - limited. 

The additional diamagnetic phase between bulk T c and 
normal state detected in temperature sweeps is marked 
by open squares. This phase transition is clearest in mag- 
netic fields between 0.2 T and 0.6 T. A linear fitting over 
these squares extrapolates down to almost at H c2 (T c ) and 
up to 1 T where the diamagnetic phase was not seen. It 
is well known that surface superconductivity boundary 
is given by H c3 (T) = 1.695F c2 (T). 15 This curve fits well 
the diamagnetic feature observed in our experiments, so 
it makes sense to assign it to surface superconductivity. 

Lastly, the strongest pinning region where theoretical 
critical current is maximum is marked by open triangles 
with a linear fitting line shown in dashed dots. It lies 
deep inside the superconducting state and is probably 
due to the matching effect. For comparison, we include 
the data for H c2 (green stars) and a crossover from vor- 
tex liquid to vortex glass phase (violet solid line), both 
estimated from the magnetization measurements. 14 The 
upper critical field is in a good agreement with our data. 
The liquid-glass crossover line, on the other hand, does 
not extrapolate to the line of the strong pinning found 
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FIG. 8. H-T phase diagram. Open and solid circles rep- 
resent bulk superconducting transition, H C 2(T), determined 
by temperature- and field-sweepings, respectively. Red solid 
line is fitting over experimental H C 2(T) using the HW theory. 
Open squares indicate the temperatures where diamagnetic 
response was detected above bulk T c . Red dashed line rep- 
resents H C 3 — 1.695H C 2(T). Open triangles represent the 
strongest pinning. Black dashed dots are a linear line fitting 
for the strongest pinning. For comparison, we include the 
data for H C 2 (green stars) and a crossover from vortex liquid 
to vortex glass phase (violet solid line) from Ref. 14. 



in our measurements (although there is no overlapping 
temperature interval). This indicates that the maximum 
theoretical critical current line is not an extension of the 
irreversibility line and represents a different crossover in 
the pinning mechanism reflected in the non-monotonic 
change in the shape of the pinning potential. 



IV. SUMMARY 

Temperature-and magnetic field-dependent penetra- 
tion depth was measured in single crystals of SrPd2Ge2 
superconductor by using a tunnel diode resonator tech- 
nique. For H = 0, the London penetration depth sat- 
urates at low temperatures indicating a fully gapped 
superconductivity in SrPd2Ge 2 . The calculated super- 
fluid density is best described by a single - gap s-wave 
superconductor, perhaps slightly on the stronger cou- 
pling side. The H c2 (T) was measured down to 0.02T C 
and could be well described by the HW theory with 
H c2 (0) = 0.4817 T and is clearly limited by the orbital 
depairing. In the magnetic fields between 0.2 T and 0.6 
T, another diamagnetic phase other than bulk supercon- 
ductivity is clearly seen, and this phase line is consistent 
with the surface superconductivity bound by the third 
critical field, H c ^ — 1.695iJ C 2- The Campbell penetration 
depth and the theoretical critical current density both ex- 
hibit non-monototic behavior with the strongest pinning 
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at intermediate temperatures, probably due to matching 
effects between the temperature - dependent coherence 
length and the relevant pinning landscape lengthscale. 

ACKNOWLEDGMENTS 

RP, MAT and HK thank V. Kogan for insightful com- 
ments. HK and NHS thank Jcehoon Kim for useful dis- 



cussions. The work at Ames was supported by the U.S. 
Department of Energy, Office of Basic Energy Sciences, 
Division of Materials Sciences and Engineering under 
contract No. DE-AC02-07CH11358. The work at GIST 
was supported by the Ministry of Education, Science and 
Technology of the Republic of Korea (2011-0028736). 



* hyunsoo@iastate.edu 
t nakhcon@gmail.com 

* chobk@gist.ac.kr 

§ tanatar@ameslab.gov 

' Corresponding author: prozorov@ameslab.gov 

1 H. Fujii and A. Sato, Journal of Alloys and Compounds 
487, 198 (2009). 

2 N. H. Sung, J.-S. Rhyee, and B. K. Cho, Phys. Rev. B 83, 
094511 (2011). 

3 E. Helfand and N. R. Werthamer, Phys. Rev. 147, 288 
(1966). 

4 N. H. Sung, C. J. Roh, B. Y. Kang, and B. K. Cho, Journal 
of Applied Physics 111, 07E117 (2012). 

5 T. K. Kim, A. N. Yaresko, V. B. Zabolotnyy, A. A. 
Kordyuk, D. V. Evtushinsky, N. H. Sung, B. K. Cho, 
T. Samuely, P. Szabo, J. G. Rodrigo, J. T. Park, D. S. 
Inosov, P. Samuely, B. Biichner, and S. V. Borisenko, 
Phys. Rev. B 85, 014520 (2012). 

6 T. Samuely, P. Szabo, Z. Pribulovo, N. H. Sung, B. K. Cho, 
T. Klein, V. Cambel, J. G. Rodrigo, and P. Samuely, Su- 
perconductor Science and Technology 26, 015010 (2013). 

7 K. Hashimoto, A. Serafrn, S. Tonegawa, R. Katsumata, 
R. Okazaki, T. Saito, H. Fukazawa, Y. Kohori, K. Ki- 
hou, C. H. Lee, A. Iyo, H. Eisaki, H. Ikeda, Y. Matsuda, 
A. Carrington, and T. Shibauchi, Phys. Rev. B 82, 014526 
(2010). 

8 J.-P. Reid, M. A. Tanatar, A. Juneau-Fecteau, R. T. 
Gordon, S. R. de Cotret, N. Doiron-Leyraud, T. Saito, 
H. Fukazawa, Y. Kohori, K. Kihou, C. H. Lee, A. Iyo, 
H. Eisaki, R. Prozorov, and L. Taillefer, Phys. Rev. Lett. 
109, 087001 (2012). 

9 N. Kurita, F. Ronning, Y. Tokiwa, E. D. Bauer, A. Subedi, 
D. J. Singh, J. D. Thompson, and R. Movshovich, Phys. 
Rev. Lett. 102, 147004 (2009). 

10 N. Kurita, F. Ronning, C. F. Miclea, E. D. Bauer, 
K. Gofryk, J. D. Thompson, and R. Movshovich, Phys. 
Rev. B 83, 094527 (2011). 

11 R. Prozorov and R. W. Giannetta, Superconductor Science 
and Technology 19, R41 (2006). 

12 R. Prozorov and V. G. Kogan, Reports on Progress in 
Physics 74, 124505 (2011). 

13 E. H. Brandt, Reports on Progress in Physics 58, 1465 
(1995). 



14 N. H. Sung, Y. J. Jo, and B. K. Cho, Superconductor 
Science and Technology 25, 075002 (2012). 

15 M. Tinkham, Introduction to Superconductivity: Second 
Edition, 2nd ed. (Dover Publications, 2004). 

16 R. Prozorov, R. W. Giannetta, A. Carrington, and F. M. 
Araujo-Moreira, Phys. Rev. B 62, 115 (2000). 

17 J. D. Fletcher, A. Carrington, O. J. Taylor, S. M. Kazakov, 
and J. Karpinski, Phys. Rev. Lett. 95, 097005 (2005). 

18 J. D. Fletcher, A. Carrington, P. Diener, P. Rodiere, J. P. 
Brison, R. Prozorov, T. Olheiser, and R. W. Giannetta, 
Phys. Rev. Lett. 98, 057003 (2007). 

19 R. T. Gordon, M. D. Vannette, C. Martin, Y. Nakajima, 
T. Tamegai, and R. Prozorov, Phys. Rev. B 78, 024514 
(2008). 

20 H. Kim, M. A. Tanatar, Y. J. Song, Y. S. Kwon, and 
R. Prozorov, Phys. Rev. B 83, 100502 (2011). 

21 E. M. Lifshitz, L. D. Landau, and L. P. Pitaevskii, Elec- 
trodynamics of Continuous Media, 2nd ed. (Butterworth- 
Heinemann, 1984). 

22 F. Ronning, N. Kurita, E. D. Bauer, B. L. Scott, T. Park, 
T. Klimczuk, R. Movshovich, and J. D. Thompson, Jour- 
nal of Physics: Condensed Matter 20, 342203 (2008). 

23 V. G. Kogan, C. Martin, and R. Prozorov, Phys. Rev. B 
80, 014507 (2009). 

24 I. Bonalde, B. D. Yanoff, M. B. Salamon, and E. E. M. 
Chia, Phys. Rev. B 67, 012506 (2003). 

25 J. R. Cooper, Phys. Rev. B 54, R3753 (1996). 

26 R. Prozorov, R. W. Giannetta, N. Kameda, T. Tamegai, 
J. A. Schlueter, and P. Fournicr, Phys. Rev. B 67, 184501 
(2003). 

27 H. Kim et al, unpublished (2012). 

28 L. L. Zhao, S. Lausberg, H. Kim, M. A. Tanatar, 
M. Brando, R. Prozorov, and E. Morosan, Phys. Rev. 
B 85, 214526 (2012). 

29 A. M. Campbell, Journal of Physics C: Solid State Physics 
2, 1492 (1969). 

30 A. M. Campbell, Journal of Physics C: Solid State Physics 
4, 3186 (1971). 

31 V. G. Kogan, R. Prozorov, S. L. Bud'ko, P. C. Canfield, 
J. R. Thompson, J. Karpinski, N. D. Zhigadlo, and P. Mi- 
ranovic, Phys. Rev. B 74, 184521 (2006). 

32 L. P. Gor'kov, Soviet Phys. JETP 10, 593 (1960). 



